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ABSTRACT 

' Context. Several newly discovered very-high-energy (VHE; E > lOOGeV) y-ray sources in the Galaxy are thought to be associated with energetic 
, pulsars. Among them, middle-aged (> lO** yr) systems exhibit large centre-filled VHE nebulae, offset from the pulsar position, which result from 
the complex relationship between the pulsar wind and the surrounding medium, and reflect the past evolution of the pulsar. 

Aims. Imaging Atmospheric Cherenkov Telescopes (lACTs) have been successful in revealing extended emission from these sources in the VHE 
' regime. Together with radio and X-ray observations, this observational window allows one to probe the energetics and magnetic field inside these 
large-scale nebulae. 

Metiiods. H.E.S.S., with its large field of view, angular resolution of < 0.1° and unprecedented sensitivity, has been used to discover a large 
population of such VHE sources. In this paper, the H.E.S.S. data from the continuation of the Galactic Plane Survey (-80° < { < 60°, 1^1 < 3°), 
together with the existing multi-wavelength observations, are used. 

Results. A new VHE y-ray source was discovered at R.A. (J2000) = 13*'56'"00\ Dec. (J2000) = -64°30'00" with a 2' statistical error in each 
coordinate, namely HESS J1356-645. The source is extended, with an intrinsic Gaussian width of (0.20 ± 0.02)°. Its integrated energy flux 
between 1 and lOTeV of 8 x 10"'" erg cm"^ s"' represents ~ 1 1% of the Crab Nebula flux in the same energy band. The energy spectrum between 
1 and 20TeV is well described by a power law dN/dE oc E"'^ with photon index F = 2.2 ± 0.2„„, ± 0.2„.j. The inspection of archival radio images 
at three frequencies and the analysis of X-ray data from i?OSAr/PSPC and XMM-NewtonfMOS reveal the presence of faint non-thermal diffuse 
emission coincident with HESS J 1356-645. 

Conclusions. HESS J1356-645 is most likely associated with the young and energetic pulsar PSR J1357-6429 (d = 2.4 kpc, r,. = 7.3 kyr and 
E = 3.1 X 10^^* erg s"'), located at a projected distance of ~ 5 pc from the centroid of the VHE emission. HESS J1356-645 and its radio and 
X-ray counterparts would thus represent the nebula resulting from the past history of the PSR J 1357-6429 wind. In a simple one-zone model, 
constraints on the magnetic field strength in the nebula are obtained from the flux of the faint and extended X-ray emission detected with ROSAT and 
XMM-Newton. Fenni-hAT upper limits in the high-energy (HE; 0.1-100 GeV) domain are also used to constrain the parent electron spectrum. 
From the low magnetic field value inferred from this approach (~ 3-4 /jG), HESS J1356-645 is thought to share many similarities with other 
known y-ray emitting nebulae, such as Vela X, as it exhibits a large-scale nebula seen in radio. X-rays and VHE gamma-rays. 

Key words. Surveys - Gamma rays: general - ISM: individual objects: PSR J1357-6429 - ISM: individual objects: HESS J1356-645 1 



1. Introduction 

The survey of the Galactic Plane conducted with the 
H.E.S.S. (High Energy Stereoscopic System) experiment 
since 2004, covering essentially the wh ole inner Galaxy 
jAharonian et al. 1l2006btlChaves et al. I2008I) . has led to the dis- 
covery of about fifty new sour ces in the very-high-energy ( VHE; 
E > 100 GeV) y-ray domain (iHinton & Hofmannll2009l) . A sig- 
nificant fraction of these sources remain without any cle ar coun- 
terpart at lower (radio and X-rays) energies ( Aharoni an et al. I 
l2008ah . Pulsar wind nebulae (PWNe), which represent the 
largest population of identified Galactic VHE sources so far, 
could natural ly account for a larg e fraction of these unclassi- 
fied sources ( Ide Jager et al. II2009I) . These bubbles of relativis- 
tic particles and magnetic field are created when the ultra- 
relat ivistic pulsar wind interacts with the surrou nding medium 
(see iGaensler&Slaiig 120061: iBucciantinil 12001 for recent re- 
views). Such interaction leads to the formation of a so-called 
termination shock, which is thought to be the site of particle 
acceleration beyond hundreds of TeV. Luminous nebulae are 
thus observed across the entire electromagnetic spectrum, in the 
synchrotron emission from radio to hard X-rays, and through 
inverse Compt on processes and potentially jfi decay from p- 
p interactions jAmato et al. Il2003t iBednarek & Bartosil^ 120031: 
iHorns et al. 1 l2006h . in the VHE domain. Radio and X-ray ob- 
servations have revealed the complex morphology of the inne r 
PWN structure at the arcsecond scale dGaensler & Slanell2006h . 
Moreover, H.E.S.S. has pro ven itself to be capabl e of mea- 
suring, in HESS J 1825- 137 dAharonian et al. Il2006el) . spatially 
resolved spectra at the tens of arcmin scale. These VHE ob- 
servations, together with the on-going Fermi-LP^ (Large Area 
Telescope) observati ons in t he high-energy (HE; 0.1-100 GeV) 
gamma-ray domain dAtwood et al. 2009), permit one to probe 
the particle spectra in these sources, and provide a unique ap- 
proach to determining the aver age magnetic field strength (see 
Ide Jager & Diannati-Ataill2009l for a discussion in this regard). 

Two classes of VHE PWNe have recently emerged, based on 
their ob servational properties: yo ung systems, such as the Crab 



nebula dAharonian et al. I l2006dh. GO.9+0.1 ([Aharonian et al.1 
l2005ah . MSH 15-52 dAharonian et all l2005bh and the newlv 
discovered H.E.S.S. sources as sociated with the Crab-like pul- 
sars of G21.5-0.9 and Kes 75 dPiannati-Atiu et al. Il2008l) . and 
the middle-aged sources (i.e. with characteristic ages Tr ^ , 



10"* yr'), as exem plified by Vela X ( Aharonian et al. Il2006cl) 



large offset centre-filled VHE sources as being the ancient neb- 
ulae from the history of the pulsar wind inside its host SNR 
(de Jager & Djannati-Atai 2009). 

This paper deals with the results from the H.E.S.S. observa- 
tions and data analysis of one of these VHE PWN candidates, 
HESS J1356-645. This source lies close to PSR J1357-6429, a 
young (T(. - 7.3 kyr) and energetic (with a spin-down luminos- 
ity E = 3.1 X 10^* erg s"') 166 ms pulsar, discove red during the 
Parke s multibeam survey of the Galactic Plane dCamilo et al. I 
12004 . At a distance of 2.4 d2A kpc (with d2.4 - d/2.4kpc), as in- 
ferred from dispersion measure, PSR J 1357-6429 may be, after 
the Crab, the nearest young (t^ < lO"* yr) pulsar known. Follow- 
up X-ray observ ations with XMM-Newton and Ch andra have 
been performed dEsposito et al. Il2007t IZavlinll2007h . Using ra- 
dio ephemerides, an indica tion of a high X-ray pulsed fraction 
(pj > 50%) was found bv lZavlinld2007h . The X-ray spectrum 
of PSR J 1357-6429 is well described by a non-thermal (magne- 
tospheric, F - 1.4 ± 0.5) compone nt plus a blackbody radiation 
at 0.16 keV dEsposito et al. Il2007h . the latter being detected in 
only a few young radio pulsars. The fraction of the spin-down 
energy channeled into the observed 0.5-lOkeV flux of ~ 7.4 x 
10"^ 4 represents one of the lowest eve r observed amongst th e 
rotation-powered pulsars. As noticed bv lEsposito et al. I d2007h . 
PSR J1357-6429 shares similarities in this regard with the Vela 
pulsar and PSR B1706-44, for which thermal emission has 
also been detected. These authors found a 3 cr marginal evi- 
dence of diffuse emissive region at the 10" scale in the XMM- 
Newton/EPlC 2-4 keV band and set a 3 cr 2-lOkeV upper limit 
on a putative PW N with Chandr a data of 2.8 x 10^^' d^^ erg s"' . 
On the other hand. lZavlinI d2007h reported the presence of a faint 
tail-like emission extending northeast at a distance of ~ 2" from 
PSR J1357-6429, with a 0.5-10 keV luminosity of 2.3 x 10^' 
d^ ^ erg s Since the proper motion of the pulsar has not yet 
been measured, the nature of this elongated feature is still under 
debate. 

The H.E.S.S. observations, data analysis and the char- 
acteristics of HESS J 1356-645 are provided in section |2] 
Results from the analysis of archival radio continuum and X- 
ray {ROSAT/PSPC and XMM-Newton/MOS) data are presented 
in section |3] together with those recently obtained from Fermi- 
LAT observations and p r esente d in the companion paper of 
iLemoine-Goumard et al. 1 d20llh . A general discussion is given 
in section|4] 



HESS J1825-137 dAharonian e t al.ll2006eh. HE SS J1718-385 



and HESS J1809-193 ( Aharoni aiiet al. Il2007bh . In the former 2. H.E.S.S. Observations and results 



case, the VHE emission is usually unresolved and centered on 
the pulsar, while in the latter case, these VHE PWNe were found 
to be significantly extended and offset from the pulsar position. 
The large differences in the measured sizes of the VHE and X- 
ray emission regions, ranging from a factor of a few to ~ 100, 
can be explained by the difference in the cool ing timescales of 
particl es radiating in the two domains (see e.g. lAharonian et al. I 
l2006eh . The evolution of the Supernova Remnant (hereafter, 
SNR) blastwav e into an in homogeneous interstellar medium 
(hereafter, ISM; Blondin et al. 2001.) an d /or the high velocity of 
the pulsar dvan der Swaluw et al. I I2OO4 may both explain these 
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' Bearing in mind that the characteristic age may not reflect the true 
pulsar age, as this is valid only when the braking index n = 3 and when 
the pulsar's initial spin period is much smaller than that observed today. 



H.E.S.S. comprises four identical 12 m diameter lACTs located 
in the Khomas Highland of Namibia at a height of 1800 m above 
sea level. Sensitive to y-rays above ~ 100 GeV, the H.E.S.S. ar- 
ray commonly achieves an event-by-event angular resolution of 
< 0.1° and a relative energy resolution of ~ 15%. The obser- 
vations discussed here were first taken between April and June 
2005, as part of the on-going H.E.S.S. Galactic Plane survey 
which now covers the band -80° < i < 60° in longitude and 
\b\ < 3° in latitude. Additionally, two sets of dedicated obser- 
vations were taken in May 2006 and March 2007 in the direc- 
tion of the pulsar P SR J 1357 -6429 in so-called wobble mode 
dPaum et al. 1 ll997h . where data are taken with an alternating 
offset from the target position of typically + 0.7°. Events within 
3° of PSR J1357-6429 and recorded when at least three of the 
four telescopes were operational were considered. After apply- 
ing quality selection to discard data affected by unstable weather 
conditions or problems related to the hardware, the resulting 
dataset was analyzed using the standard H.E.S.S. survey anal- 
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ysis scheme (lAharonian et al. I2006bh . A radius of the on-source 
region Ocut of 0.22°, a ring background region with a mean ra- 
dius of 0.8° and the hard cuts were used. These cuts include a 
minimum requirement of 200 photo electrons per shower im- 
age, paramet rized using the Hillas moment-analysis technique 
(lHillaslll996h . for y-ray selection. The zenith angles of the ob- 
servations range from 39° to 54° (with a mean value of 45°) 
inferring a mean energy threshold of 800 GeV. The acceptance- 
corrected live time is ~ 10 h at the position of the VHE emission. 
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Fig. 1. Image of the VHE y-ray excess centered on 
HESS J1356-645, smoothed with a Gaussian of width 0.13°. 
The linear colour scale is in units of excess counts per arcmin-. 
The black contours correspond to the 5, 7 and 9 cr levels for an 
integration radius of 0.22°. The black cross indicates the best 
fit position of the source centroid together with its statistical 
error The intrinsic Gaussian source size is shown by the white 
dashed circle. The dotted circle indicates the region of spectral 
extraction (see Figure |2]). The position of PSR J1357-6429 is 
marked with a green triangle. The white solid circle represents 
the 68 % containment radius of the resulting point spread 
function. 



The H.E.S.S. excess count image, smoothed with a Gaussian 
profile of width 0. 13° in order to mitigate statistical fluctuations, 
is shown in Figure[T] An excess of VHE emission is found with 
a pre-trials peak significance of 10 cr, i.e. 8.5 cr after taking int o 
account the number of trials (see e.g. lAharonian et al. 1 l2006bl) . 
The centroid position and intrinsic Gaussian width of the source 
were determined by fitting the uncorrected excess map with a 
symmetric Gaussian convolved by the nominal H.E.S.S. point 
spread function (PSF) for this dataset (with a 68% containment 
radius of ~ 0.08°). This led to a best-fit position of R. A. (J2000) 
= 13^'56'"00^ Dec. (J2000) = -64°30'00" with statistical er- 
ror of 2' for each coordinate, and an intrinsic width of 0.20°+ 
0.02° Stat, as shown by the black cross and the dashed circle in 
Figure [U respectively. Therefore, HESS J 1356-645 is clearly 
an extended source, well described by a Gaussian profile at the 
level of available statistics. 

The spectral analysis was performed on observations at less 
than 2° from the best-fit position of HESS J 1356-645, to avoid 
any systematical effect arising from larger uncertainties in the 
reconstructed energy for showers which are reconstructed far 



off-axis. The source spectrum was determined within a circu- 
lar region of 0.34° radius, as shown by the dotted circle in 
Figure [1] This region represents an ~80% source enclosure, 
which is a compromise between the optimal signal to noise ra- 
tio and the independence to the source morphology. Using the 
same cuts as for the imaging analysis, the background is esti- 
mated using the reflected- region technique, where background 
events are selected from circular off-source regions of same an- 
gular size and offset from the ob servation position as the on- 
source region dBerge et al. I l2007h . The resulting spectrum of 
HESS J1356-645 shown in Figure|2]is well described between 1 
and 20 Te V by a power-law of the form dN/dE = No (E/1 Te V)-^ 
with photon index F = 2.2 + 0.25,„, + 0.25y5 and flux normaliza- 
tion at 1 TeV No = (2.7 + 0.9,,^, + 0.4,,.^) x'lO"'- cm"- s"' TeV"' 
(X^jdof - 1.8/4, with an associated p-value of 0.77). Adding an 
exponential cutoff to the power-law does not improve the fit. A 
lower limit on the cutoff energy of 3.5 TeV at the 95% confi- 
dence level was derived. The 1-10 TeV integrated energy flux of 
8 X 10~'^ erg cm'"^ s~' re presents ~ 1 1 % of the Crab Nebula flux 
(lAharonian et al. ir2006dh . in the same energy band. 
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Fig. 2. Differential energy spectrum of HESS J1356-645, ex- 
tracted within the dotted circular region shown in Figure [T] 
Events with energies between 0.8 and 30 TeV were binned in 
order to reach a significance of at least 2 cr per resulting bin. 
The data points were fitted with a power law whose best fit is 
shown with the black solid line. Also shown are the residuals in 
the bottom panel. 



3. Search for counterparts 



3.1. General considerations 



Catalogues of Gala ctic SNRs (iGreenI l2009l) HII regions 



dPaladini et al. I 



20031). s tar-forming complexes (iRussei il 120031) 



and WR stars ( Ivan der Hucht 2001) do not list any potential 



VHE emitters in the vicinity of HESS J1356-645. However, at 
a projected distance of ~ 7' {i.e. 5 c/2.4 pc) from the centroid 
position of HESS J1356-645, PSR J1357-6429 appears to be 
the most plausible counterpart, within the context of VHE PWN 
candidates described in sectionlTj given that its spin-down lumi- 
nosity E = 3. 1 X 10"^^ erg s ' would imply a conver sion efficiency 



of Li_ioTev/E = (0.9-2.8) x 10 



-3 d^ 
"2.4- 



Moreover, ICamilo et al. I 
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Fig. 3. Radio images (in square root scale) centered on HESS J1356-645. The colour scale is such that the blue-to-red transition 
occurs at about the respective measured rms (see Table [T]i. In each case, the beamsize of the instrument is depicted by the white 
circle in the upper right corner. The position of PSR J 1357-6429 is marked with a green cross and the white contours represent the 
levels of the H.E.S.S. correlated significance, as shown in Figure[T] The black circle denotes the region within whic h flux densities 
were calculated. Left: Molonglo Observatory Synthesis Telescope (MOST) image at 843 MHz from the MGPS-2 (iMurphv et al. I 
l2007h . The grey circle denotes the position of the catalogued point-like source MGPS J135623-643323. Middle: Parkes image at 
2.4 GHz (Duncan et al. 1995|). The three grey ellipses mark the positions of the SNR c andidates G309. 8-2.6 (in the centre), G310.6- 
2.0 and G 3 10.5-3.5 (only partly sh own in the image) found bv lDuncan et al. I (Il997h . Right: Parkes-MIT-NRAO (PMN) image at 
4.85 GHz (iGriffith & Wrighli[T99l. The gr ey ellipses mark the positions of the two catalogued extended sources, PMN J 1355-6433 



and PMN .T1355-6435 (IWright et al. 111991 . 



(120041) have shown that PSR J 1357-6429 is located near the 
radio SNR candidate G309.8-2.6 (see section [X2] ). previously 
discovered and catalogued as an elon gated area of enhanced ra- 
dio emission bv lDuncan et al. I (Il997h . from the Parkes 2.4 G Hz 
survey of the Southern Galactic Plane (iDuncan et al. |[T995l see 
section [3T2I1. 



-1/3 



The distance estimate of PSR J 1357 -6429 would result in 
a VHE intrinsic source size of ~ 9 d2A pc. The observed offset 
between PSR J 1357 -6429 and the centroid of the VHE emis- 
sion could be accoun ted for as a result of either a high proper 
motion of the pulsar (Ivan der Swaluw et al. 112004. or a den sity 
gradient in the surrounding medium ( iBlondin et al.1l2001h . In 
the former case, the required transverse velocity is ~ 650 d2A 
3 km s"' , where 77.3 is the true age of the pulsar relative to its 
characteristic age of 7.3 kyr Such a velocity is re latively large 
in com parison with the mean value reported by iHobbs et al. I 
(I2OO5I) (with a mean 2D speed for pulsars with < 3 Myr 
of 307 + 47 km s"') from a sample of 233 pulsar proper mo- 
tions, but similar to what has been measured for a few PSRs. 
In the latter case, the location of PSR J1357-6429, ~ 105 doA 
pc from the Galactic Plane, i.e. beyond the thin disc known 
to be highly inho mogeneous, would argue against the scenario 
of IBlondin et al. I ( |2001,) , where the off'set PWN would be ex- 
plained by the asymmetry in the SNR geometry due to the inho- 
mogeneities in the ISM. However, the lack of available molec- 
ular line data in the direction of PSR J 1357-6429 prevents one 
from drawing firm conclusions in this regard. 

Regardless of the physical explanation for the observed off- 
set between PSR J 1357 -6429 and the VHE centroid, the sys- 
tem should have already experienced the reverse shock interac- 
tion phase, where the reverse shock has collided with the en- 
tire shock surface bounding the pulsar wind. After this phase, an 
offset PWN from the initial energetic stage of the pulsar wind 
is formed. HESS J 1356-645 would then correspond to such a 
nebula radiating in t he VHE domain. Acco r ding t o the hydro- 
dynamical model of Ivan der Swaluw et al. I (l2004l) . this phase 



occurs on a timescale of < lO"* E^j (Mpj/Mo)^^^ i^o y"^' 
where E51 is the total mechanical energy of the SN explosion 
in units of 10^' erg, M^.^- is the mass in the SN ejecta, and no is 
the ambient number density in units of cm"^. Even in the case 
of a core-collapse event (M(,y =10 M©, with E51 and no equal 
to unity), this timescale is smaller than the characteristic age of 
PSR J 1357-6429, which renders this scenario plausible. 



3.2. Radio continuum observations 

As discussed in section [T] PWNe are commonly detected in 
X-rays and radio. Interestingly enough, the above-mentioned 
SNR candidate G309.8-2.6 turns out to be located at only 
~ 0.1° away from HESS J1356-645. Archival radio im- 
ages from the Molon glo Gala ctic Plane Survey (MGPS-2) at 
843 MHz (Murphv et al. l2007h . from the Parkes 2.4GHz Survey 
(lDuncanet al.lll995h and from the Parkes-MIT -NRAO (PMN) 
survey at 4.85 GHz dGriflith & Wrightl Il993h have been in- 
spected in order to shed light on the nature of this SNR candi- 
date. The main characteristics of these three surveys are reported 
in Table □ 

Images from these three radio surveys and centered on 
HESS J1356-645 are shown in Figure [3] with the position 
of PSR J 1357-6429 denoted by the green cross and the 
H.E.S.S. significance contours from Figure [T] overlaid in white. 
Extended emission is seen in each image, partly coincident with 
HESS J1356-645. On the left image, the grey circle marks the 
position of the point-like source MGPS J 13562 3-643323, with 
a flux density of 18.2 ± 2.8 m Jy, as measured bv lMurphv et al. I 
(l2007l) . lDurican et al.l d 19971) have catalogued a list of SNR can- 
didates from the Parkes 2.4GHz Survey. The selection was based 
on the source morphology, on the absence (or faint level) of co- 
incident thermal emission by inspecting the IRAS images at 60 
yum, and on the level of polarization. Three of them, namely 
G309. 8-2.6 (at the centre of the image, 15' x 35' in size), 
G3 10.6-2.0 (on the north-western region, 45' in diameter) and 
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Table 1. Radio continuum observations towards 
HESS J1356-645. For each dataset, the centre frequency, 
the angular resolution, the noise measured over a region of 
1° in size centered on the VHE source and the nominal noise 
(in parentheses) are provided. Flux densities and errors have 
been measured within the black circle shown in Figure [3] Note 
that the flux density measured at 2.4 GHz with Parkes has been 
subtracted from the underlying emission of G3 10.5-3.5 (see 
text). 





MGPS-2 


Parkes 


PMN 


Frequency (GHz) 
rms (mJy beam"') 
PSF (FWHM, arcmin) 
Flux density & error (Jy) 


0.843 
1.8(1.6) 
0.8 
0.53 ± 0.04 


2.4 
90 (17) 

10.4 
1.5 ±0.5 


4.85 
14 (7) 
4.9 
0.54 ± 0.05 



G3 10.5-3.5 (2.7° x 3.5° in diameter), with respective flux den- 
sities of 3.9 + 1.0, 11 +3 and 19 + 5 Jy, are shown in grey 
in Figure [3] (middle) . Based on the polarized intensity images, 
iDuncan et al. suggested that the first two SNR candi- 

dates may be unrelated to the larger G3 10.5-3.5, but there is no 
distance estimate for any of these three sources. On the right 
image, the two grey ellipses denote the positions of the two ex- 
tended sources PMN J 1355-6433 and PMN J 1355-6435, with 
respectiv e flux densities of 15 8 + 11 and 140 ± lOmJy, as cata- 
logued bv lWright et al. I(ll994h . The black circle shown in Figure 
[3] represents the region of spectral extraction (with a radius of 
0.18°) and was chosen to enclose most of the diffuse radio emis- 
sion seen in the MGPS-2 and PMN images. It should be noted 
that the MOST telescope does n ot detect structures o n angular 
scales larger than 20-30 arcmin dMurphv et al. Il2007h . and that 
such large sour ces are also partially s uppressed during the PMN 
data reduction (ICondon et al. |[T993l) . This would explain why 
the larger SNR candidates G3 10.6-2.0 and G3 10.5-3.5 are not 
detected in the MG PS-2 and PMN image s, while present in the 
image processed bv lDuncan et al. 1 (ll995h at 2.4GHz. 

In the MGPS-2 image^, the flux density of the point-like 
source MGPS J 135623-643323 was removed, after summing all 
the pixels inside the black circle and correcting for the beam. 
With an average brightness (~ 3 mJy beam"') of only two times 
the local rms noise, the flux density of this extended emission at 
843 MHz (see Table [1) must be considered with caution. 

In the Parkes 2.4 GHz image^ G309.8-2.6 lies above the 
broad and faint shell-like emission of the very large SNR candi- 
date G3 10.5-3.5, which might b e responsible for muc h of the dif- 
fuse emission, as mentioned bv lDuncan et al. 1 (ll997h . Therefore 
its flux density as well as its morphology have to be taken into 
consideration. Assuming that the three SNR candidates are unre- 
lated to each other, the underlying contribution of G3 10.5-3.5 to 
the final flux density within the region of interest was estimated 
by modelling G3 10.5-3.5 as a uniform elliptical shell of 19 Jy, 
with e xternal sizes of 2.7° x 3.5° as given by iDuncan et al. I 
(Il997h and internal sizes half as large. This underlying emission 
amounts to ~ 0.48 + 0.13 Jy and was subtracted from the total 
flux density (corrected for the beam). 

In the PMN image at 4.85 GHz"*, the two catalogued sources 
have not been considered as such in the calculation, given that 
some complex diffuse objects like the one shown in Figure [3] 
(right) may have been considered as multiple sources by the 
automatical procedure, which was optimized to primarily de- 

^ [http : //www . astrop. physic s .usyd. edu. au/mosaics/| 

^ pittp : //ww w ■ atnf ■ csiro . a u/research/su rveys/2 . 4Gh_Southern/| 

"* [ftp ://ftp ■ atnf . csiro . au/pub/data/pinn/surveys/| 



tect point sources (IWright et al. 1 19941) . Instead, individual fluxes 
have been simply summed within the black circle and the flux 
density was corrected for the beam. 

The flux densities and respective errors at these three fre- 
quencies are reported in Table [T] Fitting these data points with a 
power-law gives a normalization at 1 GHz of 0.54 + 0.04 Jy and 
a slope a = 0.01 + 0.07 (where Sv cc y"), with x'/dof = 1.14/1. 
As seen in Table [1] the 2.4 GHz flux density is ~ 2 cr above the 
baseline formed by the two data points at 0.843 and 4.85 GHz. 
However, as presented above, estimating the contribution of the 
large and faint shell-like emission from G3 10.5-3.5 is a diffi- 
cult task, and the flux density measured at 2.4 GHz has certainly 
to be taken as an upper limit. For instance, the tail-like struc- 
ture north- east to the central br ight emission, originally con- 
sidered by IDuncan et al. I ( Il997h as part of G309.8-2.6, could 
simply be a region of enhanced shell emission within G310.5- 
3.5. In any case, the 2.4 GHz measurement only influences the 
normalization of the radio spectrum. A power-law fit on the 
two other data points, with a normalization at 1 GHz fixed to 
0.54 Jy, gives a consistent spectral index of ff = 0.01 ± 0.06. 
Changing the size of the region defined for calculating the flux 
densities does not change the results significantly. Such a flat 
radio spectrum suggests G309.8-2.6 to be a PWN rather than 
a shell-typ e SNR given that typica l indices for PWNe are -0.3 
< a < (iGaensler & Sland 1200 6^. while those for shell-type 
SNRs are usually close to -0.5 (Green 2009). The lack of ev- 
idence for a shell-type morphology in the MGPS-2 and PMN 
images strengthens this scenario. 

3.3. X-ray observations 

VHE PWNe may also be detected in the X-ray domain, 
as exemplified by Vela X, for which the VHE emission 
matche s quite well the extended structure seen by R OSAT and 
ASCA dAharonian et al. I l2006d: iHorns et al. I l2007l) . Such an 
extended X-ray emission has also been m easured recently 
in th e direction of HESS J 1 825- 1 37 (lUchivama et al. 
'2008!) and HESS J1809-193 (iKargaltsev & Pavlov 
,2007i) . XMM-Newton and Chandra observations towards 
HESS J1356-645 h ave been perforr ned and analyzed by 
lEsposito et al. I (l2007h and IZavUnI (l2007l) . but no clear evidence 
of a PWN has been found yet, with the exception of a faint 
tail-like emission reported by the latter group. 

3.3.1. ROSAT 

ROSAT/PSPC archival data^, with an average angular res- 
olution of 45"(FWHM, according to the Calibration Memo 
CAL/ROS/93-015), have been inspected in order to search for 
diffuse X-ray emission towards HESS J 1356-645. Figure |4] 
shows the ROSAT/PSPC image (smoothed with a Gaussian with 
cr - 5') of the vignetting-corrected and background-subtracted 
count rate in the 1-2.5 keV band, centered on HESS J1356-645. 
The observed extended and faint emission cannot be explained 
by point-like sources since only two very faint sources, flagged 
as true detections, were found within 0.2° from the centroid of 
the exten ded emission , in th e ROSAT All-Sky Survey Source 
Catalogs (IVoges etal.lll999l l200flh : IRXS J135615.9-642757 
is marked with a very low reliability, while IRXS J135605.5- 
642902 has a faint count rate of (4.81 ± 1.33) x 10"^ s"'. In 
order to calculate the flux and significance of this extended X- 
ray emission, a source region and four background regions were 

^ |http : //www . xray . mpe . mpg . de/cgi-bin/rosat/data-browser| 
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Fig. 4. ROSAT/PSPC image of the vignetting-corrected and 
background-subtracted count rate in the 1-2.5 keV band. The 
image has been smoothed with a Gaussian of width cr - 5'. 
The white contours denote the levels of the H.E.S.S. corre- 
lated significance, as in Figure [T] The position of the pulsar 
PSR J 1357-6429 is marked with a green triangle. The yellow 
dashed circles show the source and background regions used to 
estimate the significance of the extended X-ray emission. 



defined, as depicted by the yellow dashed circles of 0.2° in radius 
in Figure |4] 

After summing counts within the source (Csrc) and back- 
ground (Cbkg) regions, the excess count rate was calculated 
as Csrc -ax Cbkg, where a represents the ratio of the ex- 
posure times in the source and background regions. Since 
the ROSAT/PSPC exposure map of the region centered on 
HESS J 1356-645 is quite uniform (with a mean value of ~ 
0.45 ks and a dispersion of about 3%), a is simply 1/4. The as- 
sociated significance was computed according to the prescrip- 
tion of |Li & Ma ( 1983). The excess count rate of the X-ray dif- 
fuse emission is found to be 0.142 s"' at the 6 cr confidence 
level. Varying the sizes of the source region between 0.15 and 
0.25° does not change the results significantly. A 2D symmetri- 
cal Gaussian fit was applied to the ROSAT/PSPC image to quan- 
tify the significance of the observed source extent. The best-fit 
Gaussian width is found to be 0.14° + 0.04° (90% confidence 
level), much larger than the PSPC angular resolution. This ex- 
tended emission does not appear in the images at lower ener- 
gies (i.e. in the 0. 11 -0.4 and 0.4-1 keV bands), most likely be- 
cause of the relatively large hydrogen column density. Therefore, 
its flux has been calculated through the HEASARC software 
WebPIMMS''. With N h = (2-7) x lO^' cm"-, as measured by 
lEsposito et aP (l2007h towards PSR J 1357-6429, and assuming 
a power-law with spectral index between 2 and 2.4 (as measured 
by H.E.S.S., see also sectiongl, the ROSAT/PSPC count rate of ~ 
0.142 s"' translates into an unabsorbed 1-2.5 keV flux of 3.6^Q g 
X 10"'~ erg cm"^ s"'. 



3.3.2. XMM-Newton 

Due to the low level of statistics in the ROSAT data, no spectral 
information on the nebula could be derived. In order to constrain 
the spectral properties of the extended X-ray emission seen by 



PSPC, a 85 ks archival XMM-Newton observation (PI : G. Pavlov, 
Obsid: 0603280101), centered on PSR J1357-6429, was ana- 
lyzed using SAS v9.0.0. To remove the proton flare contamina- 
tion during the observation, a histogram of the 10-12 keV count 
rates of the two MOS cameras was built. A Gaussian fit was 
then performed on this 10-12 keV count rate distribution in or- 
der to remove time intervals where the count rates were beyond 
3 cr from the mean value (Pratt & Arnaud 2002). The remaining 
exposure time after flare screening is 55 ks. Note that, for this 
observation, the PN camera was operated in timing mode and is 
therefore not used in the analysis. The instrumental background 
was derived from the compilation of filter wheel closed obser- 
vations ^, in the same detector area as the source, and renor- 
malized in the 10-12 keV band over the whole field of view. 
Events from the list of the 42 point sources (PSR J 1357 -6429 in- 
cluded), automatically detected in the observation by the XMM- 
Newton Survey Science Center (for a total 0.5-6 keV flux of 1.4 
X 10"'' erg cm"^ s"'), were removed in order to focus on the 
properties of the diffuse X-ray emission. 
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Fig. 5. XMM-Newton/MOS image of the 0.5-6 keV flux, adap- 
tively smoothed to a signal-to-noise ratio of 15 (using asmooth 
from the SAS software) with all the detected point sources, ex- 
cept PSR J 1357 -6429, removed. Small-scale structures close to 
PSR J 1357 -6429 are better revealed in the inset image with a 
colour scale ranging from to 8 x 10"^ ph cm"^ s"' arcmin"^. 
The solid white contours denote the levels of the H.E.S.S. cor- 
related significance, while the black contours correspond to the 
ROSAT/PSPC smoothed excess count rates. The yellow dashed 
circle shows the source region used for spectral extraction, as in 
Figure |4] 



The vignetting-corrected and background-subtracted com- 
bined image of the MOS 1 and MOS 2 cameras in the 0.5- 
6keV band is shown in Figure|5] An extended emission is clearly 
seen around the position of PSR J 1357-6429, matching the ex- 
tended emission seen with ROSAT/PSPC. Small-scale structures 
are also seen close to PSR J 1357 -6429, and are discussed in 
more details, using Chandra observatio ns, in the companion pa- 
per of iLemoine-Goumard et al. 1 (l20Tl . 



|http : //heasarc . gsf c .nasa . gov/Tools/w3pimms .html| 



' |http : //xinm2 . esac . esa . int/external/xmin_sw_cal/background/| 
£ilter_closed/ 
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The yellow dashed circle shown in Figure |5] represents the 
region of spectral extraction for the X-ray nebula, which corre- 
sponds to the ROSAT/PSPC source extraction region shown in 
Figure |4] The astrophysical background was extracted over the 
rest of the field-of-view. The latter spectrum is well described by 
an absorbed APEC** model {kT = 0.21 keV, norm = 4.4 x 10"^^ 
cm-5) + power-law (F = 1.95, M) = 7.9 x 10-"* cm^^ s"' keV"' 
at 1 keV) model, with =711 for 624 dof. The hydrogen col- 
umn densit y was fixed for both components to the best-fit value 
derived by iLemoine-Goumard et al. I (1201 ih from Chandra ob- 
servations of PSR J1357-6429 A^h = (3.9 ± 0.4) x 10^' cm"". 
When let free in our model, the best-fit column density is A^h - 
(2.9 + 0.9) X 10^' cm"^ (errors are quoted at the 90% confidence 
level), similar to the value derived on the pulsar 

These two best-fit models have been renormalized to the area 
of the source region and were used as fixed components in the fit 
of the spectrum of the source region (see Figure|6l). The spectrum 
of the X-ray nebula (in red) is described by an absorbed power- 
law with the column density fixed to the value derived on the 
pulsar, F = 1.82 + 0.04 and M) = (2.65 + 0.07) x lO^^ cm"" s"' 
keV ' at 1 keV, with x' = 683 for 631 dof. The unabsorbed 1- 
2.5 keV flux then amounts to (4.25 ±0.20) x 10"'^ erg cm"- s"', 
which is in agreement with the ROSATfY'SY'C flux estimated in 
section [J. 3. 11 within the respective uncertainties. 

In order to estimate the level of systematics arising from such 
a large X-ray emission (~ 0.2° in radius, of the same order as 
the XMM-Newton field of view), several source positions around 
PSR J 1357-6429 have been tested. The resulting systematic un- 
certainties on the nebula's best-fit parameters are 0. 1 on the pho- 
ton index and 10% on the normalization. 




Energy (keV) 

Fig. 6. XMM-NewtonfMOS 0.5-10 keV specti'um of the total 
emission (in black) in the source region shown in Figure [5] The 
background spectrum (in blue) is best described by an absorbed 
power-law + APEC model, while the source spectrum (in red) 
is fitted with an absorbed power-law. A rebinning was applied 
in order to reach 10 cr/bin. Residuals are shown in the bottom 
panel. 



3.4. Fermi-LAT observations 

After nearly three years of observation, the LAT instrument on- 
board Fermi has been successful in identifying the HE counter- 
part of several PWNe powered by energetic PSRs, all of them 
being also detected in t he VHE domain (see ISlane et aL ]|20l3 
lAckermann et al. 11201 ll and references therein). The high spin- 
down flux at Earth of PSR J1357-6429 (^ 5.4 x 10^^ d,^ erg 
s~' kpc~^) is in the range of the previously LAT-detected P WNe 
(lAckermann et al. ||201 Ih . lLemoine-Goumard et al. I ( I20TTI) have 
recently discovered -y-ray pulsations from PSR J1357-6429 us- 
ing more than two years of Ferm/-LAT data, together with ra- 
dio rotational ephemerides obtained with the Parkes radiotele- 
scope. Upper limits on the emission from the associated PWN 
HESS J1356-645 in the 100 MeV - 100 GeV range have also 
been reported, by assuming a Gaussian shape with cr -0.2°, as 
measured with H.E.S.S. These 95% CL (confidence level) upper 
limits, shown in Figure |2l were derived using off-pulse data be- 
low 4.5 GeV, and from the whole dataset above 4.5 GeV, where 
no pulsed emission is detected. 

4. Discussion 

It was shown in the previous section that an extended radio 
structure, originally catalogued as a SNR candidate, and an ex- 
tended X-ray emission lie close to the newly discovered VHE 
source HESS J1356-645. Such radio-X-ray-VHE association 
is usually the observational feature of shell-type and plerionic 

* Astrophysical plasma emission code, see 

|http://hea-www.harvard.edu/APEC| 



SNRs. Regarding the morphology alone, the current VHE data 
do not allow us to firmly distinguish between these two scenar- 
ios. Despite the high confidence level of the H.E.S.S. detection, 
the VHE emission could only be characterized as a 0.2° symmet- 
ric Gaussian. The extended emission seen at 843 MHz appears as 
a faint blob of ~ 0.18° in radius with an average flux density of 
only twice as high as the rms noise, while that seen at 2.4 GHz, 
catalogued as G309.8-2.6, lies within a complex diffuse emis- 
sion, potentially associated with another SNR candidate. Only 
the radio emission observed at 4.85 GHz is safely detected. The 
ROSAT/PSPC extended emission in the 1-2.5 keV band roughly 
exhibits the same size as those in the radio and VHE domains, 
but given its significance of 6 cr, it could not be well charac- 
terized neither morphologically nor spectrally. However, thanks 
to a deep XMM-Newton observation towards PSR J 1357-6429, 
this faint and extended X-ray emission has been clearly detected. 
Therefore, the presence of a nearby energetic and young pulsar 
PSR J 1357-6429, the spectral index of the extended structure 
measured in radio and the X-ray-VHE association all point to- 
wards the scenario of a PWN associated with the past evolu- 
tion of PSR J 1357-6429. In the following, it is then assumed 
that the radio and X-ray sources are really the counterparts of 
HESS J 1356-645, and the reliability of such a scenario based 
on qualitative arguments will be discussed, keeping in mind that 
further multi-wavelength observations are needed in order to 
trigger more detailed calculations. 

In a one-zone leptonic model, a distribution of accelerated 
electrons cools radiatively by means of synchrotron (SYN) and 
inverse-Compton (IC) channels in a single volume. The IC scat- 
tering is assumed to take place in the Thomson regime. Within 
this emission volume, a constant and entangled magnetic field 
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Fig. 7. Spectral Energy Distribution of HESS J1356-645. The blue and red lines represent the SYN and IC emissions respectively, 
derived from a parent electron spectrum in the form of a power-law with an exponential cutoff (see text). The multi-wavelength data 
are presented in sections |2]and|3] The magnetic field strength is set to 3.5 jxG. 



and a homogeneous seed photon distribution from the Cosmic 
Microwave Background Radiation (hereafter, CMB) only are 
considered. For a power-law distribution of accelerated electrons 
Ky^'', the general equation relating the SYN energy flux {^syn) 
produced by electrons with Lorentz factors between yi^sYN and 
72.8 YN and the IC energy fl ux (F/c) radiated between y\jc and 
T2./C can be expressed as in lAharonian et al. 1 (l2009h : 



SYN 

Fic 



Ub 



CMB 



3-p _ 3-p \ 



(1) 



where Ug (~ 2.5 B^^ eV cm ^, where B_5 is the magnetic field 
value in units of 10 /^G) and Ucmb (= 0.26 eV cm""*) are the 
energy densities of the magnetic field and of the CMB photon 
field, respectively. It is further assumed that electrons emit syn- 
chrotron and inverse-Compton photons at the following charac- 
teristic energies (the s o-called (5-functional approximation, see 
lAharonian et al. Ill997h : 



ysYN - lAx\(fEyL,...,B-'l^ 



■"SYN.keV" -S 

yjc - 3.6x\(fE]l\ 



electrons, and for F = 2 {i.e. p = 3), equation [2l reduces to t he 
standard formula: VsynI'Pic - 10 B^^ (lAharonian et al. Ill997l) . 

Within this simple leptonic model, the VHE (H.E.S.S.) and 
X-ray {ROSATI?S?C and XMM-Newton) measurements then 
permit an exploration of the acceptable range of B-field val- 
ues in HESS J1356-645. Since ROSATjVSVC does not provide 
any spectral index, F is simply assumed to be that measured by 
H.E.S.S. TvHE - 2.2 ± 0.2. The X-ray flux measured in the 1- 
2.5 keV energy range is provided in section 13.3.11 As for the 
XMM-Newton results, the photon index F^ - 1.82 + 0.11^ is 
marginally consistent with Tvhe at the 2 cr confidence level, so 
that F is supposed to lie between 1.8 and 2. Accounting for all 
the uncertainties in the VHE and X-ray measurements, equation 
pleads to a B-field of 4.5 (2.5-8.5) yuG from 7?05Ar/PSPC and 
3.5 (2.5-4.5) juG fmmXMM-Newton. 

The extended structure seen in the X-ray and VHE y-ray 
domains would then represent the ancient PWN from the early 
phases of PSR J 1357-6429, where electrons which had been 
accelerated in the past simply age nowadays through radiative 
losses. Thus, a maximal energy is obtained by equating the PSR 
age {-2tc I (n-I), with n the braking index) to the total radia - 
tion lifetime (see equation 6 in lde Jager & Dian nati-Ataifl 2009l) . 
The associated maximal energies of IC and SYN photons then 
read: 



where "EsynmcV and YLjcjeV are the SYN and IC photon energies 
in units of keV and TeV, respectively. Substituting the above 
relations into equation [1] leads to the general expression of the 
magnetic field: 



B 



G(F)x 



SYN 



- E^-r ) 



2-r 

l.SYN.keV 



^LSYN.keV' J 



/) 



1/r 



(2) 



where G(F) ^ (0.1 x 15^"-)'^^, and F = (p + l)/2 is the common 
photon index (dN/dE oc E '^) of SYN and IC spectra. In case the 
SYN and IC fluxes both originate from the same population of 



105 yr 



0.033 



n - 1 



2 + l.KBISnG)^ 



TeV 



^IC,max 

1 TeV 



(3) 



(4) 



For a B-field in the nebula of 3.5 /iG and a PSR age of 
7.3 kyr (with n = 3), the IC and SYN photon maximal ener- 
gies would then lie at ~ 17 TeV and ~ 0.5 keV, respectively. 



' Statistical and systematical errors on the spectral index and normal- 
ization, provided in section [J.3.2l have been summed quadratically. 
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The latter value seems too low to account for the X-ray spec- 
trum as measured with XMM-Newton, well fitted with a power- 
law up to ~ 7 keV. However, it should be noted that the rel- 
evant timescale r in such E,„av estimates is rather the time 
elapsed since the end of the PWN crushin g phase, which typi- 
cally ends a few kyrs after th e pulsar birth (iBlondin et al. Il200lt 
Ivan der Swaluw et al. 1 12004 . i.e. at the time the ancient PWN 
such as HESS J1356-645 is formed. Therefore, for a B-field 
of 3.5 ;uG, equations |3] and |4] give Eic.max ~ 200 (t/2.5 kyr)"- 
TeV and EsYN.max ~ 5 (t/2.5 kyr)"^ keV, compatible with the 
XMM-Newton and H.E.S.S. measurements. Such a maximal en- 
ergy in the VHE domain, at 13.8 ± 2.3 TeV, has been de- 
tect ed in Vela X dAharonian et al. l2006d) . This measurement 
led Ide Jager & Diannati-Atail (l2009l) to constrain the magnetic 
field in the nebula to be as low as ~ 3 yuG, comparable to what 
is found here for HESS J1356-645. The VHE spectrum of the 
latter is well fitted with a power-law between 1 and 20 TeV, and 
the 95 % CL lower limit on a cutoff energy of 3.5 TeV derived in 
section|2]is compatible with these estimates. 

Within this one-zone leptonic framework, the exact expres- 
sions of SYN and IC emissions from a parent electron spectrum 
have been calculated (Blumenthal & Gould 1970) and qualita- 
tively fitted to the multi-wavelength data for a B-field of 3.5 /^G, 
as shown in Figure [T] Besides the CMB, the Galactic interstel- 
lar radiation field at the location of HESS J 1356-645, used to 
ca lculate the IC em i ssion, was derived from the latest estimates 
of iPorter & Stron j (l2005h '°. It comprises the Galactic infrared 
(from dust, at T ~ 35 K and 350 K) and optical (from stars, at T 
~ 4600 K) emission, with energy densities of 0.66 and 0.94 eV 
cm"^, respectively. In addition to the radio. X-ray and VHE y- 
ray data, the Fermi-LAT upper limits discussed in section 13.41 
dLemoine-Goumard et al. 1 1201 lb , are also shown. The parent 
electron spectrum is in the form dN^/dEg oc Eg^P exp(-Ee/Ecut) 
for Eg e [E^in, ^max], with spectral index p = 2.5, E^in - 13 GeV 
(and Emax set to 1 PeV), Ecut = 350 TeV, and a total lepton energy 
Etot = 5 X 10"*^ dj 4 erg. Note that a low energy cutoff of 13 GeV 
in the electron spectrum is in the ra nge of minimum particle 
energies considered in s everal works (iKennel & Coronirilll984t 
lAckermann et al. 11201 ih . The total particle energy must be com- 
pared to the rotational kinetic energy of PSR J 1357-6429 since 
birth, Esti = (Qq - Q^)//2, where the spin periods P,- = 27r/0,- (the 
indices and t denote the initial and current values, respectively) 
and / is the moment of inertia of the neutron star (taken to be 1 .4 
X 10"^^ g cm-). With P, = 166 ms, a lower limit on the rotational 
kinetic energy E^,; > 5 x 10**^ erg is derived for any intial period 
Pq < 70 ms. Therefore, this spin-down energy implies a reason- 
able conversion efficiency t] - Etot/E„/ <0.l 4. 

It should be noted that these calculations implicitely assume 
that the observed radio. X-ray and VHE morphologies indeed 
probe the same emissive region, whereas in Figures |4] and |5] 
there seems to be a slight offset between the ROSAT/XMM- 
Newton and H.E.S.S. extended emissions. X-rays coming from 
a more compact region closer to PSR J1357-6429, i.e. where 
the magnetic field could be higher. Only through detailed calcu- 
lations which take into account the spatial and temporal evolu- 
tion of both particles and m agnetic field in such PWN system 
(see e.g. iGelfand et al. Il20()9l) . could these key parameters be ef- 
ficiently constrained. 



5. Conclusion 

A new VHE source, namely HESS J1356-645, has been dis- 
covered during the continuation of the H.E.S.S. Galactic Plane 
Survey. This extended source lies close to PSR J 1357-6429, a 
recently discovered young (r^ = 7.3 kyr), nearby (d - 2.4 kpc) 
and energetic pulsar (E = 3. 1 x 10^^ erg s '). Archival radio and 
X-ray data have revealed an extended structure, though faint, 
coincident with the VHE emission. Given the faintness of the 
emission seen in these energy domains, the poor information 
on its morphology prevents one from distinguishing between 
shell-type and plerionic SNR as the origin of HESS J1356-645. 
However, (1) the presence of an energetic pulsar, (2) the centre- 
filled morphology with a flat spectral index measured in radio, 
and (3) the existence of a diffuse X-ray emission point towards 
the scenario of an evolved and offset VHE PWN. 

Besides comparable ages, HESS J1356-645 shares several 
similarities with the Vela PWN (whose pulsar PSR B0833-45 
has a Tc = 11.3 kyr and Hes at a distance of 290 pc), regarding 
both the pulsar (very low Lx/E efficiencies and presence of ther- 
mal X-ray emission) and the PWN (similar ratios of compact/X- 
ray to diffuse/radio PWN sizes) properti es. The 2" ta il-like X- 
ray emission close to PSR J 1357-6429 (IZavlinll2007h could be 
the pulsar jet, which is observed in the Vela pulsar at the sub- 
arcmin scales. The radio counterpart of HESS J 1356-645, ~ 
0.35° in size, would thus correspond to the bright 2° x 3° ra- 
dio structure (the so-called Vela X), extending southwest of the 
Vela pulsar In terms of the broad-band non-thermal emission, 
the so-called Vela cocoon, a 45'-long collimated X-ray filament 
extending so uthward of PSR B0833- 45, has been detected in 
VHE y-rays dAharonian et al. 1 120063) . while the whole Vela X 
radio str ucture was recently de tected in the GeV domain with 
AGILE dPelUzzoni et al. 1 1201 Oh and Fermi-ENY dAbdo et al. I 
l20T0bl) . These measurements seem to point towards the sce- 
na rio of two distinct lepton spectral com ponents, first suggested 



by |de Jager (2001) and la ter studied by Ide Jager et arr(12008h 



|http : //galprop . Stanford. edu/resources . php?option=data| 



and iLaMassa et alii d2008l) . However, follow-up H.E.S.S. ob 
servations have recently allowed for the detection of signifi- 
cant VHE emission beyond the cocoon, up to 1.2° from the 
VHE barycentre, wi th a spectrum similar to that of the cocoon 
dPubois et al. II2009I) . In the case of HESS J 1356-645, at the es- 
timated distance of PSR J1357-6429, the putative "cocoon-like" 
structure would appear as nearly point-like (~ 0.09° long) for 
H.E.S.S. However, such structure cannot be the dominant VHE 
component as HESS J 1356-645 is well fitted simply with a 
Gaussian width of 0.2° (see section|2li. HESS J 1356-645 would 
then represent the VHE counterpart to the "Vela X-like" ra- 
dio structure (see section 13.2b . detected also in X-rays with 
ROSATI^SVC and XMM-Newton (see section O- Thus, given 
the observational limitations of the current radio and VHE y- 
ray data, the broad-band spectrum of HESS J 1356-645 can be 
explained by a single lepton population with reasonable param- 
eters, as presented in section |4] 

From qualitative arguments, it was shown that the magnetic 
field within HESS JI356-645 must be quite low (~ 3-4 pG), 
of the same order as what has been im plied in other PWNe 
{e.g. Vela PWN, see lde Jager et al. l2008h . under similar assump- 
tions as investigated here. Follow-up observations in radio will 
help to constrain the pulsar braking index and proper motion, 
and confirm the nature of the extended emission. The combined 
GeV- TeV data from the continuation of Fermi observations and 
future VHE observations with the upcoming H.E.S.S. phase II 
will better constrain the spectral shape of the IC emission, in 
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order to probe the magnetic field strength and the energetics of 
HESS J1356-645. 
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